Time–frequency analysis for music signals

Time–frequency analysis for music signals is one of the applications of time–frequency analysis. Musical sound can be more complicated than human vocal sound, occupying a wider band of frequency. Music signals are time-varying signals; while the classic Fourier transform is not sufficient to analyze them, time–frequency analysis is an efficient tool for such use. Time–frequency analysis is extended from the classic Fourier approach. Short-time Fourier transform (STFT), Gabor transform (GT) and Wigner distribution function (WDF) are famous time–frequency methods, useful for analyzing music signals such as notes played on a piano, a flute or a guitar.

1 Knowledge about music signal
2 Short-time Fourier transform
- 2.1 Continuous STFT
- 2.2 Discrete STFT
3 STFT example
4 Spectrogram
5 Wigner distribution function
- 5.1 Formula
6 See also
7 Sources

Knowledge about music signal

Music is a type of sound that has some stable frequencies in a time period. Music can be produced by several methods. For example, the sound of a piano is produced by striking strings, and the sound of a violin is produced by bowing. All musical sounds have their fundamental frequency and overtones. Fundamental frequency is the lowest frequency in harmonic series. In a periodic signal, the fundamental frequency is the inverse of the period length. Overtones are integer multiples of the fundamental frequency.

Table. 1 the fundamental frequency and overtone
Frequency	Order
f = 440 Hz	N = 1	Fundamental frequency	1st harmonic
f = 880 Hz	N = 2	1st overtone	2nd harmonic
f = 1320 Hz	N = 3	2nd overtone	3rd harmonic
f = 1760 Hz	N = 4	3rd overtone	4th harmonic

In musical theory, pitch represents the perceived fundamental frequency of a sound. However the actual fundamental frequency may differ from the perceived fundamental frequency because of overtones.

Short-time Fourier transform

Continuous STFT

Short-time Fourier transform is a basic type of time–frequency analysis. If there is a continue signal x(t), we can compute the short-time Fourier transform by

$\mathbf{STFT} \left \{ x(t) \right \} \equiv X(t, f) = \int_{-\infty}^{\infty} x(\tau) w(t-\tau) e^{-j 2 \pi \tau} \, d \tau$

where w(t) is a window function. When the w(t) is a rectangular function, the transform is called Rec-STFT. When the w(t) is a Gaussian function, the transform is called Gabor transform.

Discrete STFT

However, normally the musical signal we have is not a continuous signal. It is sampled in a sampling frequency. Therefore, we can’t use the formula to compute the Rec-short-time Fourier transform. We change the original form to

$X(n \, \Delta t,m \, \Delta f) = \sum_{p=n-Q}^{n%2BQ} x(p \, \Delta t) e^{-j 2 \pi p m \, \Delta t \, \Delta f} \, \Delta t$

Let $t = n \, \Delta t$ , $f = m \, \Delta f$ , $\tau = p \, \Delta t$ and $B = Q \, \Delta t$ . There are some constraints of discrete short-time Fourier transform:

$\Delta t \, \Delta f = \frac{1}{N},$ where N is an integer.
$N \ge 2Q%2B1$
$\Delta < \frac{1}{2f_\max}$ , where $f_\max$ is the highest frequency in the signal.

STFT example

Fig.1 shows the waveform of a piano music audio file with 44100 Hz sampling frequency. And Fig.2 shows the result of short-time Fourior transform (we use Gabor transform here) of the audio file. We can see from the time–frequency plot, from t = 0 to 0.5 second, there is a chord with three notes, and the chord changed at t = 0.5, and then changed again at t = 1. The fundamental frequency of each note in each chord is show in the time–frequency plot.

Spectrogram

Figure 3 shows the spectrogram of the audio file shows in Figure 1. Spectrogram is the square of STFT, time-varying spectral representation. The spectrogram of a signal s(t) can be estimated by computing the squared magnitude of the STFT of the signal s(t), as shown below:

$\mathbf {spectrogram} (t,f) = \left| \mathbf{STFT} (t,f) \right|^2$

Although the spectrogram is profoundly useful, it still has one drawback. It displays frequencies on a uniform scale. However, musical scales are based on a logarithmic scale for frequencies. Therefore, we should describe the frequency in logarithmic scale related to human hearing.

Wigner distribution function

The Wigner distribution function can also be used to analyze music signal. The advantage of Wigner distribution function is the high clarity. However, it needs high calculation and has cross-term problem, so it's more suitable to analyze signal without more than one frequency at the same time.

Formula

The Wigner distribution function $W_x(t,f)$ is:

$\mathbf W_x(t,f) = \int_ {-\infty}^\infty x(t%2B\tau/2)x^*(t-\tau/2) e^{-j2\pi\tau\,f} \,d \tau,$

where x(t) is the signal, and x*(t) is the conjugate of the signal.

Sources

Joan Serra, Emilia Gomez, Perfecto Herrera, and Xavier Serra, "Chroma Binary Similarity and Local Alignment Applied to Cover Song Identification," August, 2008
William J. Pielemeier, Gregory H. Wakefield, and Mary H. Simoni, "Time-frequency Analysis of Musical Signals," September,1996
Jeremy F. Alm and James S. Walker, "Time-Frequency Analysis of Musical Instruments," 2002
Monika Dorfler, "What Time-Frequency Analysis Can Do To Music Signals," April,2004
EnShuo Tsau, Namgook Cho and C.-C. Jay Kuo, "Fundamental Frequency Estimation For Music Signals with Modified Hilbert–Huang transform" IEEE International Conference on Multimedia and Expo, 2009.